Ferromagnetism in the Hubbard model on fcc-type lattices 
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The Hubbard model on fcc-type lattices is studied in the dynamical mean-field theory of infinite 
spatial dimensions. At intermediate interaction strength finite temperature Quantum Monte Carlo 
calculations yield a second order phase transition to a highly polarized, metallic ferromagnetic state. 
The Curie temperatures are calculated as a function of electronic density and interaction strength. 
A necessary condition for ferromagnetism is a density of state with large spectral weight near one 
of the band edges. 
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More than 30 years ago, the Hubbard model was in- 
troduced to describe band magnetism in transition met- 
als, in particular the ferromagnets Fe, Co and Ni [0-1). 
However, for a long time it appeared to be a generic 
model rather for onii-ferromagnetism than for ferromag- 
netism. At half filling (one electron per site) on bipartite 
lattices antiferromagnetism emerges in both weak and 
strong coupling perturbational approaches and, in par- 
ticular, antiferromagnetism is tractable by renormaliza- 
tion group methods Q. By contrast, ferromagnetism is 
a non-trivial strong coupling phenomenon which cannot 
be investigated by any standard perturbation theory. In 
fact, our knowledge on the possibility of ferromagnetism 
in the Hubbard model is still very limited. Only in a 
few special cases was the Hubbard model proven to have 
ferromagnetic order. The first rigorous result showing 
a fully polarized ground state, the theorem by Nagaoka 
1^ , is only valid for a single hole added to the half filled 
band in the limit of infinite on-site repulsion. In the 
limit of low electronic density a saturated ferromagnetic 
ground state has been found recently in the one dimen- 
sional Hubbard model on zigzag chains . Further, the 
existence of a ground state with net polarization has been 
proven for the half filled band case on bipartite lattices 
with asymmetry in the number of sites per sublattice , 
and in 'flat -band' systems [^. 

Due to the development of new numerical algorithms 
and powerful computers the problem of ferromagnetism 
in the Hubbard model has recently become accessible to 
numerical investigations of finite systems, at least in re- 
duced dimensions (d = 1,2). In the presence of a next 
nearest neighbor hopping, polarized ground states were 
found very recently in d = 1 |]l0| and on square lattices 
in the case of a van Hove singularity at the Fermi energy 

The rigorous and numerical results mentioned above 
show that the stability of ferromagnetism is intimately 
linked with the structure of the underlying lattice and the 
kinetic energy (i.e. the hopping) of the electrons. This 
fact has also been observed by exact variational bounds 
for the stability of saturated ferromagnetism |p^-p^ , and 
within approximative methods |l6). Generally, lattices 
with closed loops and with frustration of the competing 



antiferromagnetism (non-bipartite lattices) are expected 
to support ferromagnetism. Non-bipartite lattices have 
an asymmetric density of states (DOS) and thus a peak 
away from the center of the band. Indeed, it has already 
been observed by the inventors of the model that a 
DOS with a peak near one of the band edges is favorable 
for ferromagnetism. Hence, the fee lattice is expected to 
be a good candidate for ferromagnetism because i) it has 
a highly asymmetric DOS and ii) antiferromagnetism is 
frustrated. 

Since numerical studies of finite systems would allow 
only very small fee lattices, we study the Hubbard model 
in the dynamical mean-field theory (DMFT) which be- 
comes exact in the limit of infinite dimensions, d = oo 
We will show that an appropriate, non-singular 
kinetic energy is able to induce metallic ferromagnetism 
in the single band Hubbard model for intermediate values 
of the interaction. 

In the limit d — s- oo, the system is reduced to a dy- 
namical single site problem ||l9(| equivalent to an Ander- 
son impurity model complemented by a self-consistency 
condition [^,^. Still, it cannot be solved analytically, 
and to avoid further approximations we employ a finite 
temperature Quantum Monte Carlo method ||2|] . Unlike 
traditional mean-field theories, the DMFT takes quan- 
tum fluctuations fully into account. Spatial fluctuations 
are neglected - an approximation which becomes exact in 
the limit of large coordination number, Z. For the three 
dimensional fee lattice we have Z = 12. The DMFT has 
proven to be a powerful and reliable tool for the study 
of three dimensional correlated Fermi systems both for 
thermodynamical as well as dynamical properties p8| , p3[ . 
Spectral properties are obtained by analytic continuation 
of the imaginary time data using the Maximum Entropy 
method The solution of the self-consistent single 

site problem provides the local self energy Eo-n = ^uii^n) 
{ll}„ ■ Matsubara frequencies) from which the local Green 
function is obtained by a Dyson equation: 
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= iix)n + ^J■~ Scrn- Here, e(k) are the single par- 
ticle energies, No{E) is the non-interacting DOS, and L 
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is the number of lattice sites. Investigations of the Hub- 
bard model on the hypercubic lattice in d = oo did not 
show indications for ferromagnetism for finite U Jl8| , p3t , 
however in a recent study at [/ = c» a ferromagnetic 
instability was found within the non-crossing approxi- 
mation 25 1 . 

A non-trivial generalization of the fee lattice in d di- 
mensions [0 (see also [|6|) is the set of all points with 
integer cubic coordinates summing up to an even integer. 
It is a non-bipartite Bravais lattice for any dimension 
d > 2. Each point has Z — 2d{d — 1) nearest neighbors 
defined by all vectors x which can be written in the form 
X = ±ei ± ej , with two different cubic unit vectors ei and 
Bj — 1, - ■■ ,d). The energy dispersion reads 



e(k) = 4t 



cos ki cos kj 



l<i<j<d 
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where we allow a next nearest neighbor hopping t' . In the 
particular case t' = t/2 we find e(k) = t e\^{k) — 3t where 
e?ic(k) is the energy dispersion of the hypercubic lattice. 
This relation leads to an inverse square-root divergency 
of the DOS at the lower band edge in any dimension. 
In the limit d = oo this divergency occurs for any t' , 
and with the proper scalingi = 1/Vz the DOS can 
explicitly be calculated as |14| : 
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was also used by Uhrig's |2^ calculation of the 
single spin-flip energy of the fully polarized state. There, 
the fcrromagnet turned out to be stable against a single 
spin-flip over a wide density regime. 

In d = 3 for < < t/2 the divergency is absent, and 
the DOS has a typical ^/E — i?;,— behavior at the lower 
edge Et = e(7r, 7r/2, 0) = -4t + 2t' . For i' = there is a 
weak, logarithmic divergency at E},. The main effect of t' 
on the DOS is to induce a broad peak with much spectral 
weight near (but not right at) the lower edge. In (^) and 
in the following the energy scale is set by the variance of 
the non-interacting DOS. In d = 3 for < t' < t/2 the 
total band width isW ~ 16i + 4<'. By fixing the variance 
of the DOS, V \/l2i^T6?2 = the band width ranges 
from 4.618 for i' = to 4.899 for t' = t/2. Comparing 
with real 3d-transition metals, e.g. Ni, one can roughly 
identify our energy scale with leV. 

We will consider the following two cases: i) strictly 
d = oo, i.e. using the DOS (||) and ii) d = 3 lattices, 
i.e. using the DOS of the the three dimensional lattice in 
(|l|) . In practice we perform finite sums over the k- vectors 
corresponding to a finite three dimensional lattice. To 
keep the finite size error at least one order of magnitude 
smaller than the statistical errors, the number of lattice 
sites (respectively k-vectors) has to be of the order of 10^ 
to 10^ depending on temperature. The computer time for 
this summation is still much smaller than for the Monte 
Carlo sampling. 



To detect a ferromagnetic instability we determine the 
temperature dependence of the uniform static suscep- 
tibility, xf, from the two particle correlation functions 
p7| . In addition, we calculate commensurate magnetic 
susceptibilities and the compressibility (charge suscepti- 
bilty) . No charge instability is observed in the parameter 
regime under consideration. 

At an intermediate interaction strength of L/ = 4 we 
find the ferromagnetic response to be strongest around 
quarter filling (n ~ 0.5). xf obeys a Curie-Weiss law 
(see Fig. 1 for case i)) and the Curie temperature Tc can 
safely be extrapolated from the zero of Xp^ to a value of 
Tc = 0.051(2) at n = 0.58 (2|]. Below the magneti- 
zation M grows rapidly, reaching more than 80% of the 
fully polarized value {Mmax = n = 0.58) at the lowest 
temperature which is only 30% below Tc. The three data 
points M (T) (Fig. 1) are consistent with a Brillouin func- 
tion with the same critical temperature of Tc = 0.05 and 
an extrapolated full polarization at T = 0. A saturated 
ground state magnetization is also consistent with the 
single spin-flip energy of the fully polarized state |M . 
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FIG. 1. Magnetization, M (diamonds), and inverse uni- 
form static susceptibility, Xf^ (circles; values multiplied by a 
factor of 4 to use the same scale) for U — 4 and n — 0.58. 
Error-bars are of the size of the symbols or smaller. Dotted 
line is a linear fit to Xf^' dashed line a fit with a Brillouin 
function to M. (Note that the circle at T = 0.05 is a data 
point, not an extrapolation.) 

The Curie temperatures are obtained from the Curie- 
Weiss plots for different values of interaction U and elec- 
tronic densities n = Nd/ L, leading to the phase diagram 
Fig. 2 for the case i). At low temperatures the system 
is ferromagnetic over a wide density regime. Tc increases 
with U and the upper critical density Uc seems to ap- 
proach n = 1. The values of function of U are 
consistent with the single spin-flip results |^ . Note that 
the Stoner criterion, UN{fj,) > 1, would give Uc > 1 for 
U > 1.5. For low densities the Curie temperature be- 
comes very small and seems to vanish close to n = 0. 
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Antiferromagnetism is not expected on the fee lattice in 
high dimensions because the difference of the numbers of 
not frustrated bonds and frustrated bonds is only of the 
order of d resulting in an effective field of the order of 
t^d(xl/d M. 
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FIG. 2. Tc vs. n phase diagram of the Hubbard model on 
the d — oo fee lattice with [7 = 2 (full circles), (7 = 4 (open 
circles), and U = 5 (squares). Lines are guides to the eye. 
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FIG. 3. Phase diagram of the Hubbard model on the d = 3 
fee lattice with t' = t/4 and U — 6. Circles: ferromagnetic, 
squares antiferromagnetic phase boundary. 

In case ii) we first studied the pure fee lattice, i.e. t' — 
0. However in this situation, and for intermediate [/—val- 
ues (U < 6) the Curie-Weiss extrapolation would lead 
to very small Curie temperatures, numerically indistin- 
guishable from Tc = 0. If we, however, introduce a next 
nearest neighbor hopping term of i' = i/4 the situation 
is again similar to the case i). The phase diagram for 
U ^ 6 (Fig- 3) shows a large ferromagnetic regime with 
an optimal density about n — 0.65. For U — A no fer- 
romagnetism was found for temperatures accessible to 
the Monte Carlo technique. In contrast to the case i) 



now there seems to exist a finite lower critical density. 
This might be due to the fact that the DOS has no di- 
vergency at the band edge but a foot. An additional 
feature is the antiferromagnetic phase very close to half 
filling. We found an instability at k = (7r,0, 0), corre- 
sponding to ferromagnetic planes with alternating orien- 
tation in direction perpendicular to the planes. In addi- 
tion there seem to be instabilities at k = (7r/2, 7r/2, 7r/2) 
and k = (7r/2, 7r/2, 0), however only at n = 1 and very 
low temperatures (T < 0.01). At n = 1 the Neel tem- 
perature is (numerically) degenerate with an instability 
at k = (7r,7r/2,0). Below half filling the latter is more 
strongly suppressed. 

The single particle spectra for spins parallel and an- 
tiparallel to the net magnetization (Fig. 4 for case i)) are 
apparently metallic since both spectra have a finite value 
at the Fermi level {w — ii). This also holds in the param- 
agnetic phase (not shown). The minority spin spectrum 
(dotted line in Fig. 4) is not only shifted to higher fre- 
quencies but has also a little foot at low energies. This 
foot contains about 50% of the spectral weight below /i 
and looses weight with increasing polarization. The mi- 
nority spin spectrum shows a pronounced upper band 
around w — /i « [/ = 4 with a developing (pseudo) gap. 
The majority spin spectrum, however, is only slightly 
affected by the interaction. Here, the weight of the up- 
per band is very small since the Pauli principle makes 
it unlikely for the majority spin-electrons to hop over 
occupied sites. 
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FIG. 4. Single particle spectrum, N{uj), for (7 = 4, 
n = 0.58, and T = 0.04. FuU (dotted) line: majority (mi- 
nority) spin spectrum. 

The exchange splitting, defined by the difference of 
the peak positions, is A « 0.43 which is significantly 
smaller than the shift A'^*""'^'' = UM w 2 in Stoner 
theory. Further, Stoner theory would give a Curie tem- 
perature T^sto"'^'- « 0.5 which is one order of magnitude 
too high. Thus, the mean-field Stoner theory cannot de- 
scribe the different energy scales correctly. Identifying 
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our energy scale with leV, the resulting Curie tempera- 
ture corresponds to Tc « 600A'. We note that this value 
and the large ratio A/Tc > 8 - which will even some- 
what increase for T — *■ - are in good agreement with 
experimental values on fee Ni {Tc = 627 K, A/T^ « 10) 
|30| . The agreement is quite surprising considering the 
neglect of band degeneracy and inter-orbital interactions 
favoring ferromagnetism (Hund's coupling). 

While the saturated ferromagnet has zero interaction 
energy {U{d) = 0), the mobility of the electrons and 
hence the kinetic energy is very poor, since electrons can- 
not hop over occupied sites. The paramagnetic state, on 
the other hand, also tries to minimize double occupan- 
cies, (d), and thus has a small kinetic energy, too. For 
example at U = 4 and n = 0.58, we find U{d) < 0.01 for 
T < 0.06. Hence the gain in interaction energy cannot 
account for the relatively high Curie temperatures. This 
raises the question: under which circumstances has the 
fully polarized ferromagnet a lower kinetic energy than 
the paramagnet? In a one-particle picture at T = the 
kinetic energy difference is: 



AE 



dEENo{E)~2 dEEN{E) (4) 



where the Fermi levels Ep^rro and Ep have to determined 
by the density n. In N{E) is the (unknown) DOS 
of the correlated paramagnet. To get an idea of the en- 
ergy difference we consider the following simple picture: 
since double occupancies are effectively suppressed by U , 
the band will consist of two parts separated by an en- 
ergy U . For a given t-electron the average number of 
nearest neighbors which are not occupied by |-electrons 
is q = 1 — ni = 1 — n/2in the paramagnet, and hence 
the effective width of the lower band will be reduced by 
that very factor q. The simplest assumption for the lower 
part of the spectrum is N{E) = No^E/q) i.e. a renormal- 
ization of the energies (Sl). For the DOS (|), AE{n) is 
negative for all densities < n < 1 and vanishes for n = 
and n = 1. The density dependence of |A_E(n)| resembles 
strongly the Tc(n) data (Fig. 2). In particular |A£'(n)| 



becomes maximal at n = 0.735 with lAE'L 



0.0747. 



For a DOS with the same shape but a divergency at the 
upper edge, AE is always positive, i.e. the ferromagnet is 
totally unstable. Apparently this simple treatment of the 
correlated paramagnet describes qualitatively the stabil- 
ity of the ferromagnet depending on the structure of the 
non-interacting DOS, and gives even the correct energy 
scale for Tc- To check this approximation one can also 
estimate AE by replacing the second integral of (^ by 
a numerical integration over the finite temperature DOS 
in the paramagnetic solution. For n = 0.58, U = A, 
and T ~ 0.04 (cf. Fig. 4, but now in the paramagnetic 
solution) we obtain AE' « —0.04(1) which is in rea- 
sonable agreement with the simple approximation above 
{AE = -0.064). 

In summary, we presented a numerical proof of metal- 
lic ferromagnetism in the single band Hubbard model in 



the dynamical mean-field theory. A density of states 
with large spectral weight near one of the band edges is 
an essential ingredient for ferromagnetism. This condi- 
tion goes far beyond the Stoner criterion, UNq{ij,) > 1, 
because it is the band renormalization of the interact- 
ing paramagnet - a manifest many body effect - which 
determines the stability of ferromagnetism. 
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